Enhanced cooperativity below the caging temperature of o-terphenyl 
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The utility of a cooperative length scale for describing the dynamics of glass-forming 
liquids is shown using literature data on o-terphenyl. Molecular dynamics and Monte Carlo 
simulations reveal a distribution of cooperative fractal events below the caging temperature 
Ta- Guided by these results, we show how to extract the size of slow regions in any glass- 
forming liquid from probe rotation/diffusion measurements, which agrees quantitatively with 
4-D NMR and grows steadily as temperature is lowered below Ta- We clarify why this length 
must also be the size of the largest cooperative events. 
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Many liquids either cannot crystallize or crystallize sufficiently slowly that they vitrify below 
their glass transition temperature T g . A fundamental understanding of glass formation is still 
lacking because it has not been firmly established whether the pronounced slowing down is simply 
kinetic in origin or there is an underlying thermodynamic character Building upon 

previous experimental results, this letter shows a natural length scale for cooperative motion in 
o-terphenyl that grows below a caging temperature, giving strong support to the thermodynamic 
viewpoint Jl|. The caging temperature Ta, is seen in many dynamics experiments on glass-forming 
liquids Jj, \A , \^ . Particularly lucid are experiments that measure translation and rotation of small 
probe molecules immersed in the liquid. Above Ta, these are coupled by the Stokes-Einstein relation 
and have the same temperature dependence. In contrast, below Ta these measurements decouple, 



with translation having a weaker temperature dependence than rotation 



A wide 



variety of ideas have been expounded to explain this observation 



12, 131 and these ideas fall 



into two categories. One group of ideas suggest that translation is dominated by a different and 



faster moment of a distribution of relaxations than is measured by rotation 



n 



The second 



12j, |l3j, and in this paper we present 



group of ideas consider both to be controlled by slow motions 
compelling evidence that this is the correct picture. 

n 

The 1965 model of Adam and Gibbs 1] suggests that there should be cooperative motion in 
glass- forming liquids. The cooperative size £, is inversely related to the configurational entropy of 
the liquid S con f. At temperatures sufficiently above T g , all molecules undergo independent local 



2 



Brownian movements without signs of cooperativity, where the relaxation time r a , and the viscosity 
7/, vary as 



r a ~ 77/T ~ exp [E/k B T] T > T A 



(1) 



where E is the high-T activation energy. As temperature is lowered, the density of the liquid 
gradually increases and Brownian motion becomes hindered, as neighboring molecules block each 
others attempts to move. This crowding leads to cooperative dynamics [X P. llX Ilfil. . active 
for all T below the caging temperature Ta- The onset of cooperativity is also accompanied by 
the observed 'caging effect' in the mean-square displacement of a particle between the ballistic 
and diffusive regimes, and a reduction in S con f, causing £ to grow. These changes result in a 
progressively stronger temperature dependence of 77 and r Q at lower temperatures. 

Since experimental attempts to identify the length scale for cooperative motion have met with 



limited success 

Q,0,Q,Q, 
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6|, the dominant evidence for this quantity is from computer simulations 

^nn 
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23L |2jJ . Simulations have the profound advantages of direct observation 



of motion and straightforward identification of both the size and shape of cooperatively rearranged 
regions. Equilibrium simulations of liquids are not yet possible near T g , but have been done down 
to 0.77a [22|. In addition to confirming the essential aspects of the Adam and Gibbs model, 
simulations have provided two novel insights. Instead of a sin gle s ize scale for cooperative motion, 
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26]. 



I here is in fact a broad distribution of si:j stales below 7 4 [15I H, Q, Q, 
The largest size in this distribution £ grows rapidly as temperature is lowered, as expected by 
Adam and Gibbs |l|. The second important observation is that the cooperatively rearranging 
regions are not the three-dimensional volumes that were initially proposed, but instead are fractal 
21, Q. The observed fractal dimension (of order 2) clearly shows that the majority 
of molecules within the volume £ 3 have not participated in the cooperative event. Consequently 
a new model for cooperative motion was proposed that accommodates these new insights |27j ]. 
All glass-forming liquids show a temperature dependence of cooperative size, with commensurate 
effects in viscosity and relaxation time, in reasonable accord with the expectations of dynamic 



scaling 



27, 



28 



£ ~ [(T - T C )/T C ] 



-3/2 



T C <T <T B 



(2) 



r a ~ £ 6 exp(E a /k B T) ~ (T — T c )- g exp(E a /k B T) 



(3) 
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T/T A 



FIG. 1: Simulation results for the length scale of cooperative motion as a function of reduced temperature. 
From molecular dynamics simulations of Lennard- Jones sphere mixtures (▲ j^l) and from Monte Carlo 
simulations of the bond fluctuation model at density 0.8 for chains of 10 monomers Inset is the 

temperature dependence of the configurational entropy (o |22() and the product of cooperative length and 
configurational entropy (□ 26] ) for chains of 10 monomers. 



r, ~ £ 6 exp(E v /k B T) ~ (T - T c )- 9 exp(E v /k B T) (4) 

The critical temperature Tq, is below T g and is the temperature at which the equilibrium extrapo- 
lated values of £, r a and rj all diverge. The material-specific activation energy apparently depends 
on whether 5 e gm e„ tal relation ( E „) or „ (E v) is nreasured. with E v ,E a , 1.2 for „o„- 
polymeric organic glasses [29(. Tb is the crossover temperature p, |30|] below which all molecules 
are caged by their neighbors (T g <Tb < Ta)- Tg is typically about 20 K below Ta and corresponds 
to the upper temperature limit where dynamic scaling (and Vogel-Fulcher) quantitatively describe 
the temperature dependence of relaxation. 

Figure H shows data for £ as obtained from molecular dynamics simulations [3] and new results 



22 



. In the MC simulations 



from Monte Carlo (MC) simulations on the bond fluctuation model 
the mobile particles are identified as those that move over a time scale of interest. We then look 
for clusters of these mobile particles and find that their sizes are a function of time for short times, 
but quickly become time independent up to r a . We only consider cluster sizes in this intermediate 
time range. Figure ^ shows the sizes of the largest mobile particle clusters as a function of T. Both 
sets of simulations show that £ is sensibly independent of temperature above Ta, but then grows 
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rapidly when temperature is lowered below Ta- The very different nature of the simulations used 
for the data in Fig. ^ suggests there may be universal aspects to the growing length scale below 
Ta- The inset of Fig. ^ shows that S con f as defined by Adam and Gibbs also changes character 
at Ta with a broad crossover between Tb \^ (below which Eqs. 01 and 0] or Vogel-Fulcher should 
describe dynamics) and a much higher temperature above which dynamics obey Eq. ^ j^lt - The 
inset also shows that the product of £,S con f is essentially temperature independent. Hence, the 
growing length scale is reciprocally related to the vanishing configurational entropy, in qualitative 
agreement with Adam and Gibbs llll. 

n 

Guided by simulations, and owing to the abrupt change in the very nature of relaxation at Ta [la 
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34 1, the caging temperature is easily identified by a variety of experiments probing liquid 



dynamics. We demonstrate this point with o-terphenyl using rotation and translational diffusion 
of molecular probes, dielectric spectroscopy, viscosity and Fabry-Perot interferometry. Molecular 
probe techniques provide a model independent measure of Ta- The length scales extracted from 
these experiments, which are in quantitative agreement with existing £(T) data from 4-D NMR 
show a strong temperature dependence only below Ta- 4-D NMR 

Q, 0, Q, Q, Q, Q 

provides the benchmark length scale of the slow regions, but unfortunately this method can only 
be used over a very limited temperature range j^ . 

Tb was found to be 286K from extrapolation of the /3-relaxation time to the a-relaxation 
(Fig. |2j A), in agreement with the literature 2^1 . The temperature dependence of r a and rj (Fig. El A) 
were fit using Eqs.|3]and0]and values for Tq, E a , and were determined to be 227 K, 120 kJ/mol 
and 140 kJ/mol respectively. In figure El A, equation ^ was fit to n using E = 22 kJ/mol. This 
activation energy was then used to fit r (with a result that extrapolates to about 10 -14 sec at 
an infinite temperature). A Ta of 313 K was then defined by the crossover between Eq. ^ and 
Eqs.E]and|U The Ta from this crossover between Eqs. 01 and 0] and Arrhenius behavior is also the 
temperature below which the Stokes-Einstein relation breaks down. This 'breakdown' below Ta is 
apparent in figure ElB where both probe and self-diffusion show a temperature dependence below 
Ta I? , 47 , 4||- The calorimetric T g of o-terphenyl was previously determined to be at 241 K 49]. 

The Stokes-Einstein relation expects the rotational relaxation time (r r ) and translational dif- 
fusion coefficient D of probe molecules are coupled so that their product is independent of tem- 
perature. Above Ta the probes diffuse a distance of order their own size in the time it takes for 
the probe to rotate. However, at Ta these two dynamics decouple, a fact which can be used to 
establish the cooperative length scale. To understand this, consider a bimodal distribution with a 
small quantity (volume fraction (j)f « 1) of fast and <p s = 1 — (j)f slow particles j^]. The average 
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B. Model Independent 
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FIG. 2: Dynamic data for o-terphenyl above and below Ta ■ Model dependent measures of T A (A): 
viscosity (o |42t l43| ~l . dielectric spectroscopy (a-relaxation (A 44]) and /3-relaxation (0 |4a])) and Fabry- 
Perot interferometry (V Q]). The dashed line is Ta- The dotted line is a fit to the /3-relaxation data 
extrapolated to Tb- Solid curves are fits to t q and r\ using Eq. Q]for T > Ta and Eqs.[3]and2]for T < Tb- 
Model independent quantities affected by Ta (B): self-diffusion (x and tetracene probe dynamics (+ 
HQ)- Solid curve is E q-EI 

rotational time is dominated by the slow molecules 

( Tr ) = (j) s T s + 4>fTf = 4> s t s (5) 
since <j) s <j)f and t s ^> tj. The diffusion coefficient is the sum of fast and slow contributions 

D = 4> S D S + 4>fD f = <M s 2 /6r s + ij> f S 2 f /&r, (6) 



where ^ = 6D s t s and £| = 6-D/T/. Since 
be used to define a length scale, 



1, the product 6L>(r r ) 



e s + 



&r./Tf 



C = 6 = V®D(T r ). 



can 



(7) 



Fig. |3lA shows the temperature dependence of this length scale for two probe molecules in o- 
terphenyl. We find that the length scale defined by Eq. quantitatively agrees with 4-D NMR, 



which is known to target the slow contribution £ s . This quantitative agreement requires 4>iCs ^ 
4>f4>s£ 2 f T s /Tf, justifying Eq.[3 

In the calculation of 6 over a wide temperature range, (r r ) was interpolated using the noted 
relationship (r r ) ~ r\jT 9| along with the viscosity data of Fig. |2J ^ from rotation/translation of 
anthracene (circles in Fig. EJ) shows a strong temperature dependence over the entire experimental 
temperature range because all T < Ta, whereas the larger tetracene probe (triangles in Fig. 03) 
gives £ that is clearly temperature dependent below T A and independent of temperature above Ta- 
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FIG. 3: The temperature dependence of the cooperative length scale in o-terphenyl (Tq = 227 K). From 
4-D NMR (■ with error bars [3]) and calculations using probe dynamics (o anthracene j^J; A tetracene 



[4JJ). Dashed lines denote T g = 241 K, T B = 286 K and T A = 313 K. The solid line is the slope of -3/2 
expected by dynamic scaling (Eq. [3J • 



Ideally a probe molecule would have the same shape, size, polarity and properties as the liquid 
matrix in which it was inserted. ^-NMR can provide measurements of translational and rotational 
se//-diffusion coefficients |2L This allows for direct measurement of £ using Eq. Although 
measurements of (r r ) exist jflj they are rare. In the case where measurements of (r r ) are lacking, 
we assume that rj/T properly describes the temperature dependence of (r r ) [?], 13] ■ This assumption 
allows for measurements of Ta that stand in agreement with the other measurements of Ta (see x 
symbols in Fig. [2J, but with the limitation of not being able to give an absolute measure of £. 

In Fig. |31 measurements of £ are plotted along with the absolute measurements of £ provided 
by 4-D NMR (£ = 2.9 ± 1 nm at 252 K 36]). While each technique measures £ differently, they 
all stand in quantitative agreement. This result strongly suggests that the vast majority of the 
molecules are "slow", <j>f ~ 10~ 5 , which makes the slow term dominate even if t s /tj = 10 3 . This 
concept ^of the slow molecules dominating near-T 9 dynamics has also been suggested by recent 



models 



12 



3- 



The cooperative length estimated from figure is 5.2 ± 1.0 nm at T g and 0.33 ± 0.06 nm at 
Ta- The van der Waals sphere radii r v ^w OI o-terphenyl was calculated from atomic radii using 
the procedures of Edward ^]. A r v dw = 0.37 nm for o-terphenyl is within the calculated range 
of £(Ta), suggesting that the magnitude of cooperative size calculated herein is reasonable. This 
result £(Ta) = r V( iw appears to be general [3]. 

The length scale for cooperative motion has been estimated from measurements of probe diffu- 
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sion and rotation. At T g + 10 K this method agrees quantitatively with 4-D NMR and also provides 
a measure of the cooperative size at higher temperatures. This length decreases as temperature 
is raised and adopts the van der Waals radius of the glass- forming liquid above Ta- The size of 
the slow regions exhibits the temperature dependence expected by dynamic scaling for the largest 
cooperatively moving (fast) regions, suggesting that there is only one important length scale in 
this problem. This makes sense because cooperative sizes need not grow any larger than the size 
of the slow regions. 

Despite the fact that the cooperative volume is fractal instead of space-filling, the essential 
features of the Adam-Gibbs model Q] are correct. The cooperative size does indeed grow rapidly 
as temperature is lowered below the caging temperature. Over a temperature range extending from 
Ta to T g , £ of o-terphenyl has been observed to increase by an order of magnitude from the van 
der Waals radius of each molecule, making the insight of Adam and Gibbs particularly noteworthy. 

We thank the National Science Foundation (DMR-0422079) for funding and S.K. Kumar for 
discussions. 
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